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Motivated by recent experimental findings, we have developed a theory of the superconducting 
state for 2D metals without inversion symmetry modeling the geometry of a surface superconducting 
layer in a field-effect-transistor or near the boundary doped by adsorbed ions. In such systems the 
two-fold spin degeneracy is lifted by spin-orbit interaction, and singlet and triplet pairings are 
mixed in the wave function of the Cooper pairs. As a result, spin magnetic susceptibility becomes 
anisotropic and Knight shift retains finite and rather high value at T = 0. 

PACS numbers: 74.20.-z, 73.20.At, 71.18.-hy, 76.60.Cq 
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In this presentation we consider the energy spectrum 
and some other properties of the superconducting state 
for a two-dimensional (2D) system without inversion 
symmetry when the two-fold degeneracy of the electron 
energy spectrum is lifted by spin-orbit (SO) interaction. 
We have been motivated by several recent experimental 
discoveries. Superconductivity (SC) with critical temper- 
ature Tc = 91 K was observed in insulating WO3 when 
its surface was doped by Na+ ions 0J^. According to 
Rcfs. I^J^, SC in the material is a surface effect which 
comes to the existence in a narrow layer with a com- 
position of Nao.osWOa. We speculate that strong near- 
surface electric fields should be present in the SC layer 
of these crystals quite similar to the 2D superconducting 
fuUerene and polyacene crystals in the field-effect tran- 
sistor geometry The two- fold degeneracy of the 
electron energy spectrum in the bulk characteristic of 
crystals with 3D inversion symmetry should be broken 
in the surface layer because of SO coupling. This effect 
is known and well documented experimentally |5|j^ in 
physics of semiconductor microstructures, and has been 
recently discovered in surface states of metallic Au 0] 
and also of (Li covered) Mo and W In what follows, 
therefore, we discuss general properties of 2D supercon- 
ductors with lifted spin degeneracy. 

Ideal boundary being a strong "defect" for the crys- 
tal's bulk, it is known that split-off 2D electron bands 
may form at the surface (Tamm levels). Besides, near 
the surface the space inversion (the 3D inversion which 
ensures two-fold degeneracy of the bands throughout the 
Brillouin zone |^]) is automatically broken, resulting in 
a new SO invariant of the group Coov (and its discrete 
sub-groups) in the free electron Hamiltonian 



H^o = a.{(T X p) 



(1) 



a unit vector normal to the surface {h = 1). Despite the 
fact that p is two-dimensional, the 3D group controls the 
symmetry of the system since n determines the prefered 
direction of the normal to the plane and spins are always 
3D. Total single electron Hamiltonian, H = Hq + Hso, di- 
agonalizes in eigenfunctions 77a (p) and eigenvalues e\ {p) 



Vxip) = -= 



-( ' 

^2 V «Aexp(i(^p) 

exip) = eoip) + el°{p), 4°(p) = -aXp , 



(2) 



dubbed in the literature as Rashba term [ pX)|pd| ]. Here p 
is a 2D quasimomentum, a are Pauli matrices, and n is 



A = ±1, and the Fermi surface splits into two sheets. SO 
interaction Hgo puts electron spins into the plane, (_L n), 
where they are aligned perpendicular to the quasimomen- 
tum p. The magnitude of the coefhcient a depends on 
the electric field, presence of adatoms at the boundary, 
atomic weight, and atomic shells involved |l2). Angle re- 
solved photoemission spectroscopy (ARPES) data on the 
s-p surface band of Au are in excellent agreement with 
Eq. (^). They result in a nearly isotropic spectrum with 
a « 5 X 10~^ eV- cm (that is somewhat larger than in 
typical semiconductor systems) and in the SO splitting 
at the Fermi level 2|q;|pf = 0.11 eV Q. Anisotropic 
d-derived bands of Mo and W show even larger SO split- 
tings up to 0.13 eV and 0.5 eV, respectively, depending 
on the surface concentration of adsorbed Li"*" ions. 

Returning back to the surface SC reported in Refs. 
for doped tungsten bronzes, we expect, therefore, that 
the SO interaction of Eq. (|l|) should play essential role 
in the surface properties both of the normal and SC 
phases. In what follows we address changes in properties 
of the BCS-typc 2D superconductivity caused by lifting 
the electron spin degeneracy in the presence of the SO 
interaction. 

For the simplicity sake we consider a two-particle in- 
teraction C/(|r — r'l) (here and below all vectors lie in the 
plane). Retardation effects such as due to the phonon 
mediated attraction are known to be of no importance 
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for a weak coupling scheme. To put it briefly, we do not 
address at all the issue of the mechanism responsible for 
pairing, and restrict ourselves by features pertaining to 
the reduced dimensionality (for the proper symmetry at 
the surface) and spin-split isotropic electron spectrum. 
The SO splitting at the Fermi level, 2apF, (here and be- 
low a > 0) can be comparable to, or exceed, the SC char- 
acteristic scale, the critical temperature Tc- However, it 
is assumed to be small compared to the Fermi energy, 
2apF ^ ef, and p+ — p- ^ pp, where p\ — p± are 
Fermi momenta for the two SO split sheets of the Fermi 
surface. Therefore, near the Fermi energy the electron 
dispersion is taken linear: 

S.\{p)^vf{p-P\), p\ = [1 + a\/vF)pF- (3) 

The total electron Hamiltonian H = Hq + Hgo + i?int 
includes the interaction Hamiltonian 

X^ivp pp'q 

X (P)«^(-P - q)a>>(-p' - q)ap(p')- (4) 

Here ax{p) are Fermi electronic operators, and 
U\fj,up{p, p', q) is the interaction potential written in the 
representation of the spinors rjx (p) 

Ux^..pip , p',q) = C/(|p-p'|) 

X (??A(p),?7p(p'))('?M(-p-q)''?"^(-p'-q))' (5) 

Z7(|p|) is the Fourier image of the interaction potential, 
and scalar products of spinors are equal to 

{dr + ex{p))g\ip,T - t ) - \ ^^Uxxt,t,{p,p' 

MP' 



ivxip), Vx'ip')) = Hl + XX'eM-i{ipp-ipp,)]}. (6) 

In what follows we employ the method of thermody- 
namic (Matsubara) Green functions Onset of the 
SC state manifests itself as the appearance, side by side 
with the regular Green functions [Q 

gxip, t-t')^ -{T{ax{p, r)a+(p, r')}), (7) 

also of the non-zero anomalous (Gor'kov) averages 

/a(p, t-t') = X{T{~ax{p, t)5a(-P, r')}), (8) 



/+(p,r-r') = A(rK(p,r)5+(-p,r')}), 



(9) 



where operators aA(p,T) are in the Heisenberg represen- 
tation. Functions /a(p) = /a(p,0-|-) that obey the rela- 
tion 

/a(p) = /a^(-p) (10) 

are components of the wave functions of Cooper pairs at 
each spectrum branch. Fermionic commutation relations 
result in 



/a(p) 



-/a(-p) 



(11) 



From equations of motion for the operators aA(p,T), 
c'rfiA(p,''') = [H,ax{p,T)], equations of motion for 
.gA(p,T — t') and fxiP,T — t') follow in the usual way. 
They include 4-fermion T-products, the standard decou- 



pling [^|T5| of which results in 
0)U-p')f+{^p,r-r')^-Sir-r') 



(12) 



{dr - exip))f^{-p, T - r') + \Y,^^U^t.xxip', P, 0)/+(p').9a(p, r - r') = 0. 



hp' 



(13) 



For the Cooper pairs with the zero total momen- 
tum forming the condensate, q = in J7aa^^(p, p', q)- 
Next simplification in the potential comes from the fact 
that it enters in Eqs. (|lj;|l^) only in the combination 
with the functions /^(p), hence, because of Eq. (11), 
Ux\pp{p,p' ,0) should be antisymmetrisized. Simple al- 
gebra shows 



Uxt,{p,p') = ^[Uxxpp{p,p' ,0) - Uxx,^,^{p, -p' ,0)] 

« iA/iC/(0)exp[i(<^p, -<^p)]. (14) 

Having in mind s-pairing we neglected the momentum 
dependence of the potential in the right hand side of 
Eq. (p^). Making Fourier transformation of Eqs. (^2|jl^) 
to the Matsubara frequencies icun, we find the final form 
of Gor'kov equations 



[iujn - exip)]gx{p,i^n) + A(p)/+(-p,w„) = 1, (15) 

A+(p).gA(p,c^n) + + £a(p)]/a+(-P,w«) = 0, (16) 
where the "gap" A(p) is 

A(p) = exp(-i</jp)Ao, Ao = ^U{0) ^ exp{iipp)fxip), 

Ap 

A+(p) = A*(p). (17) 

Before turn to the discussion of the nature of the order 
parameter in the new SC state, let us make one prelim- 
inary comment. The "gap function" A(p) depends on 
p through its phase. This dependence is inherent in the 
non-perturbative character of the spinor basis functions 
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of Eq. (g) after spin degeneracy is lifted. It cannot be 
fully eliminated, however, it can be changed by a differ- 
ent choice of phase factors in Eq. (^. Green functions 
diagonal in branch indices iQ are 

/+(-p,L^„)- A+(p)/i^„A(p), 

A.A(p)=c.f, + e2(p) + |A|2. (18) 

The self-consistency condition for |A(p)| = A(T) is 

-^tu{0)J2d-1{p) = i. (19) 

In the approximation 2apF ^ £f, left hand side of 
Eq. ( p^ does not depend on a and coincides with the 
BCS equation for a = 0. The energy spectrum found 
from the poles of g\{p, ujn) and f\{p, uj^) with the proper 
account of Eq. (0) is 



i?A(p)-±[4(p-PA)' + A2(T)]l/2, 



(20) 



i.e., it consists of two gapped branches. One sees that 
weak SO coupling makes no changes in thermodynami- 
cal properties of the new BCS like state. In particular, 
the relation A(0) = {tt/j)Tc ~ 1.76Tc holds in this case. 

Green function were found above in the representation 
of energy branches indices A. To get insight on the sym- 
metry of the new SC state we rewrite them as matrices 
G{p,LUn) and F(p,a;„) in the basis of the components 
77^ (p) of the spinors 77^ (p) 



Gapip,UJn) ^J2^"(P)9xiP,^n)VxiPy 
A 

Fc,f3ip,uj„) = ^77j(p)/A(p,w„)?7f (p), 



(21) 



and split g\{p,uJn) and fx{p,uJn) into their symmetric 
and antisymmetric in A parts 



5a = 5s + ^A^as, f\ = fs + h>^fa: 



(22) 



Eq. (I22D implies that G = Gs + Gas and F = Fs + i^s. 
Antisymmetric parts of g\ and /a (and also of G and F) 
originate because of the splitting of the energy spectrum 
by the SO interaction. It is convenient to write Gg and 
Gas (and also Fg and Fas) in terms of the Pauli matrices 
cr and the unit matrix erg • After some algebra we get 

Gs(p, W„) = 5s(p,t^n)CTo, (23) 



Gas(P,W„) = ^5as(P,'^n)(p" X cr) • n, (24) 



Fs(p,c^„) = [/s(p,c^„)e^'^>'](Tj„ (25) 



^as(p,^n) = -[fUP,^n)e"^-]<Jy{(r-p"), (26) 

where p" = p/p. One sees that phase factors e"^*^P^ drop 
out from both Fg and Fas according to Eqs. (^Q. 
Gas and Fas components of the Green functions are of 
the order of unity but only inside a narrow ring of the 
momentum space, Ap « — p_, between two Fermi 
surfaces. 

Both contributions to G are invariants of the group 
Coov, however, they are of opposite parity if the trans- 
formation p — > — p, rotation by tt around the z axis, is 
performed only in coordinate space with spin remaining 
unchanged. Both Gas and Fas depend on the azimuth of 
p'^. The factor (cr -p") is a pseudo-scalar of the 3D inver- 
sion group. We attribute to Gas a number of transport 
effects in normal conductors with lifted spin degeneracy 
which have been observed or predicted. 

Existence of non-zero F-function, Eqs. (|25| ) and (|2^), 
means the broken gauge symmetry in the SC state, 
U{1). Eq. ([ill), when expressed in terms of the func- 
tions F(j^(p, LUn), reflects permutational symmetry of the 
wave function of the two electrons of the pair, (p, a) 
and (— p,/3). If there is the 3D inversion center, one 
can make the total wave function odd either by choos- 
ing an antisymmetric spin wave function while the space 
part is symmetric (singlet pairing, S = 0), or, vice versa, 
by choosing a symmetric spin wave function while the 
space part is antisymmetric (triplet pairing) {S = 1), 
see Ref. [|6|. Therefore, F^^ cx {<Jy)ai3 for 5 = and 
Faf3 (X {(Jy{(T ■ p))q/3 for = 1. For the broken 3D 
inversion near the surface, the Cooper pair wave func- 
tion F becomes a mixture of singlet and triplet pair- 
ings. Fs(p,a;,i) is the singlet part component, while 
Fas(p,ti'ri) provides for the triplet admixture. Note that 
this term comes about even for the short range interac- 
tion U{\y - r'l) cx (5(r - r'). 

As an example of the phenomenon in which breaking 
the parity of the order parameter manifests itself, we cal- 
culate the spin susceptibility x^'^iF) of the surface SC 
state. It can be directly tested by Knight shift experi- 
ments. For systems with strong SO coupling, the total 
magnetic susceptibility x cannot be split into the orbital 
and spin parts. However, when a <C wf such a division is 
justified ||l^, and wc find x*'^ neglecting the quantization 
of the orbital motion. Following Ref. ||l^, we calculate 
the tensor x^ {hj — Hi ^) as a linear spin response to 
the magnetic field B 



sp 



= -^|Tzy(£F)X! j tr{crjG(p, W„)crjG(p,CJ„) - atF{p,UJn)cTjF+{p,UJn)}- 



(27) 
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Here the trace is taken over the spin indices, i^{£f) = 
prf^TTVp is the density of states at the Fermi level (per 
one spin orientation), and /Ub is the Bohr magneton. In 
the normal state the right hand side of Eq. ( p7| ) diverges 
because of Landau singularity Q , and to resolve it the 
summation over ujn should be performed first. Finally 



sp 

Xn 



(28) 



Remarkably, Xn found in the weak field limit, /ieS <C 
apF, when Zeeman splitting is small compared to the SO 



splitting, is isotropic and coincides exactly with the Pauli 
susceptibility. Since the difference in x'^p in the normal 
and SC states comes from the integration over the region 
of only about A(T) near sp, the appropriate integral con- 
verges. Therefore, in these terms the integration over f 
can be performed first. Non-diagonal components of x''^ 
are equal zero because of the Cqov symmetry, while for 
the diagonal components of x^p somewhat cumbersome 
calculations result in 



1-nT 



E 



A2(r) 



(29) 



xTiT)^kxIiT) + 



1 sp 
2Xn 



E 



7rrA2(T) 



K + A2(r)) 



3/2 



(30) 



It follows from Eqs. (H) and (|^) that xT{T) ^ 
xjj^(T') and x'i(O), X||'^(0) 7^ 0, i.e., spin succeptibil- 
ity is anisotropic and does not turn into zero at T = 
0. We attribute this fact to the admixture of the 
triplet state to the s-type SC ground state because 
of the broken parity. For app <C A(T), it follows 
from Eqs. (|2|,||) that xT(0) ~ l[apF / A{T)fx''S and 
X||P(0) w i[apF/A(T)]2x^. Hence, for a = the BCS 
result xTC-*) — X||'^(0) = is recovered. In the op- 
posite limit A(T) <C app, we get xTC-*) ~ Xn ^"^^ 
xjj^(O) ~ hence, X± (^) is nearly T independent 

while xt^FiT) drops twice when T changes from to 
T = 0. 



In conclusion, we have shown that spin-orbit interac- 
tion lifts the spin degeneracy for 2D (surface) supercon- 
ductor resulting in the two gapped branches in the energy 
spectrum. Thermodynamics of such a state would be al- 
most identical to a BCS like superconductor. It turns out 
that due to broken 3D inversion symmetry at the surface, 
the pair wave function is the mixture of singlet and triplet 
components. To the best of our knowledge, such a mix- 
ture has never been studied before. As the result, strong 
anisotropy appears in the spin-susceptibility tensor and 
the Knight shift. Splitting of two gapped branches in the 
momentum space depends on the strength of SO cou- 
pling. The latter can considerably exceed the SC scale 
Tc to be observed by other means. 
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